Characterizing the nonlocal correlations of particles that never interacted 
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Quantum systems that have never interacted can become nonlocally correlated through a process called en- 
tanglement swapping. To characterize nonlocality in this context, we introduce local models where quantum 
systems that are initially uncorrelated are described by uncorrelated local variables. While a pair of maximally 
entangled qubits prepared in the usual way (i.e., emitted from a common source) requires a visibility close to 
70% to violate a Bell inequality, we show that an entangled pair generated through entanglement swapping will 
already violate a Bell inequality for visibilities as low as 50% under our assumption. 



It is natural to expect that correlations between distant par- 
ticles are the result of causal influences originating in their 
common past — this is the idea behind Bell's concept of lo- 
cal causality (H]. Yet, quantum theory predicts that measure- 
ments on entangled particles will produce outcome correla- 
tions that cannot be reproduced by any theory where each sep- 
arate outcome is locally determined by variables correlated at 
the source. This nonlocal nature of entangled states can be 
revealed by the violation of Bell inequalities. 

However remarkable it is that quantum interactions can es- 
tablish such nonlocal correlations, it is even more remark- 
able that particles that never directly interacted can also be- 
come nonlocally correlated. This is possible through a process 
called entanglement swapping j2]. Starting from two indepen- 
dent pairs of entangled particles, one can measure jointly one 
particle from each pair, so that the two other particles become 
entangled, even though they have no common past history. 
The resulting pair is a genuine entangled pair in every aspect, 
and can in particular violate Bell inequalities. 

Intuitively, it seems that such entanglement swapping ex- 
periments exhibit nonlocal effects even stronger than those of 
usual Bell tests. To make this intuition concrete and to fully 
grasp the extent of nonlocality in entanglement swapping ex- 
periments, it seems appropriate to contrast them with the pre- 
dictions of local models where systems that are initially un- 
correlated are described by uncorrelated local variables. This 
is the idea that we pursue here. To precise it further consider 
the general scenario depicted below. 




A source Si sends particles to Alice and Bob, and a sepa- 
rate source S2 sends particles to Charles and Bob. All parties 
can perform measurements on their system, labeled x, y and z 
for Alice, Bob and Charles respectively, and they obtain out- 
comes denoted a, b, and c, respectively. Bob's measurement y 
might correspond to a joint measurement on the two systems 
that he receives from each source. The correlations between 
the measurement outcomes of the three parties are described 



by the joint probability distribution P(a, b, c\x, y, z). An en- 
tanglement swapping experiment is clearly a particular case 
of this scenario, where Bob's measurement corresponds to a 
Bell measurement entangling Alice's and Charles's particles. 

Under the usual assumption, the tripartite distribution 
P(a, 6, c\x, y, z) would be said to be local if it can be writ- 
ten in the factorized form 

P(a,b, c\x,y,z) 

r dXp(X) P(a\x, \)P(b\y, \)P(c\z, A) , (1) 



where the variable A with distribution p(\) describes the joint 
state of the three systems according to the local model, and 
P(a\x, A), P(b\y, A), P(c\z, A) are the local probabilities for 
each separate outcome given A. 

In our scenario, however, they are two separate sources Si 
and S2. It it thus natural to assume that the local model as- 
signs two different states Ai and A2, one to each source, and 
to consider instead of (Q]) the decomposition 

P(a,b,c\x,y,z) (2) 
= //dAidA 2 p(Ai, A 2 ) P(a\x, Xi)P{b\y, X u \ 2 )P(c\z, A 2 ). 



The local response function for Alice now depends only on Ai, 
the one of Charles only on A 2 , while the one of Bob depends 
on both Ai and A 2 . So far, the decompositions (Q~|) and (0 
are equivalent because p(Ai , A 2 ) could be different from zero 
only when Ai = A 2 = A to recover ([]]). We now introduce 
our basic assumption: since the two sources Si and S 2 are 
supposed to be independent and uncorrelated, it is natural to 
assume that this property carries over to the local model. The 
variables Ai and A 2 should therefore be independent and their 
joint distribution p(Ai, A 2 ) factorize: 



p(Ai, A 2 ) = pi(Ai)p 2 (A 2 ). 



(3) 



We refer to models satisfying this independence assumption 
as "bilocal" models, since they aim at explaining the correla- 
tions P(a, b, c\x, y, z) with two independent sources of local 
variables. 

Even though the local variables Ai and A 2 are initially in- 
dependent, once conditioned on the joint measurement re- 
sult of Bob they will bear enough correlations to reproduce 
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non-trivial correlations between Alice's and Charles's system. 
These correlations, however, are much weaker than those that 
can be established through joint measurements in quantum 
theory. We introduce below a (quadratic) Bell inequality that 
is satisfied by all bilocal correlations, but which is violated by 
quantum correlations. As we will see, our inequality simpli- 
fies the requirements for the demonstration of quantumness in 
entanglement swapping experiments. 

Restricted classes of local models with independent sources 
were considered in within the context of the detec- 

tion loophole. But apart from these exploratory works, little 
was known about how nonlocality is induced through mea- 
surements on independent quantum systems. Beyond its fun- 
damental interest, nonlocality is also known to play a key 
role in several quantum information protocols yl |g], and 
measurement-induced correlations are at the basis of quan- 
tum repeaters |Q] and measurement-based quantum computa- 
tion 18[]. One of our contribution is to introduce a theoretical 
framework to address broadly the role of nonlocality in such 
contexts. 

Before entering in the details of our results, it might be 
worth justifying further our independence assumption. It is 
strictly speaking an assumption, rather than something which 
follows logically from locality. Indeed, some events in the 
common past of the sources S\ and S2 could in principle have 
influenced, in a way that is perfectly in accord with locality, 
both Ai and A2 such that they wind up correlated, in violation 
of (|3). However, an assumption similar to (0 is actually hid- 
den in any standard Bell-type experiment. In order to derive a 
Bell-type inequality, one needs (in addition to the premise of 
local causality) an assumption having to do with the measure- 
ment settings being "freely chosen". What this means in prac- 
tice is that the measurement settings are determined by a ran- 
dom mechanism that is considered independent of the variable 
A describing the particle source. Seen from this perspective, 
the assumption that the laser sources in the quantum random 
number generators used to choose the measurement settings 
in a standard Bell experiment |9] are independent of the laser 
source generating the entangled photons is not much different 
from the assumption that the two laser sources (which may be 
of different brands, assembled in different parts of the world, 
and powered by different electrical supplies) used in an entan- 
glement swapping experiment are independent. Of course we 
cannot exclude in principle that such apparently independent 
sources are significantly correlated. But, quoting Bell, "this 
way of arranging quantum mechanical correlations would be 
even more mind-boggling than one in which causal chains go 
faster than light. Apparently separate parts of the world would 
be deeply and conspiratorially entangled" [Q]]. 

Characterization of the bilocal set. We start by giving 
a characterization of the set of bilocal correlations which is 
more handy for analytical and numerical purposes than the 
definition (0 and ([3]). First note that without loss of gener- 
ality the local response function P(a\x, Ai) of Alice can be 
taken to be deterministic, i.e., such that it assigns a unique 
measurement output a to every input x (any randomness used 



locally by Alice can always be thought of as being included 
in the shared variable Ai). In the case of a finite number 
of possible measurement inputs and outputs, there is a finite 
number of such deterministic strategies corresponding to an 
assignment of an output a x to each of Alice's N possible in- 
puts x. We thus label each of these strategies with the string 
a = ai . . . &n and denote the corresponding response func- 
tion P a (a\x). Similarly, the response functions P(b\y, Ai, A2) 
and P(c\z, A2) can also be taken deterministic. We label the 
associated strategies (3 and 7 and the corresponding response 
functions Pp(b\y) and P 7 (c|z). 

Let A^L, denote the set of pairs (Ai,A2) specifying the 
strategies a, f3, and 7 for Alice, Bob, and Charles. Defining 
q a pj = JJ ' A i2 dAid\2 p(Ai, A2), Eq. (01 can then be rewrit- 
ten as 

P(a, b, c\x,y, z) = q a[)l P a (a\x)P(3{b\y)P 7 (c\z) (4) 

with q al 3 7 > and J2 a p-y 1<*M = ^ So far we have not 
used the independence condition (0, and (0 corresponds to 
the well-known decomposition of local correlations as a con- 
vex sum of deterministic strategies, where the weights q aPl 
can be understood as the probabilities assigned by the source 
to the strategies a, (3, and 7. 

Let us now define q ai = J2p Qapy, Qa = q a [3~ ( , and 
1i = Ec/j 1a/3 7 - Using the fact that Aj t 2 /37 = (A* x A 2 ) n 
Ai 2 , as follows from (0, the independence condition (0 im- 
plies that 

q ai = q a q~(- (5) 

Conversely, any correlation P(a, b, c\x, y, z) satisfying (0 
and (0 can be written in the form (0. Indeed, since 
q ai = q a q~<, we can write q aPl = q a q^qp\ al . In- 
serting this expression in (0 and defining P an (b\y) = 
Y^pqi3\a-iP[){b\y), we then find that P(a,b,c\x,y, z) = 
Ea,-, 9a?7 P a (a\x)P a . 7 (b\y)P 7 (c\z), which is clearly of the 
form (0. We thus conclude that a tripartite correlation is bilo- 
cal if and only if it admits the decomposition (0 with the re- 
striction (0. 

The bilocal set that we have just characterized is clearly 
contained in the local set. The extremal points of the local 
set, corresponding to deterministic strategies, are also bilo- 
cal, but a mixture of deterministic strategies is not necessar- 
ily bilocal due to the non-convex constraint (0. Therefore, 
one cannot use standard Bell inequalities to distinguish one set 
from the other. As we will see below, however, correlations 
can be shown to be non-bilocal using non-linear inequalities 
(or joint sets of linear inequalities). Note that while decid- 
ing if a given set of correlations is local can be conveniently 
solved using linear programming, deciding if a correlation is 
bilocal is a quadratically constrained problem which is much 
more difficult to handle numerically. However, standard linear 
and semidefinite relaxations can be used to deal with the non- 
linear constraint (0. We describe in Holl . a linear relaxation 
of (0 which works well on many instances. 
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Application to entanglement swapping. We now illustrate 
how the bilocality constraint restricts the set of possible corre- 
lations on a simple example inspired by the standard entangle- 
ment swapping protocol. The sources S% and 5 2 send pairs of 
particles in the singlet state |\fr _ ) = (|01) - ] 10)) /y/2. Bob 
performs a Bell state measurement on the two particle he re- 
ceives, with four possible outputs b = bobi = 00, 01, 10, 11 
corresponding to the four Bell states |$ + ), |\f r+ ), l^ - ), and 
\&~), respectively. Depending on Bob's result, Alice and 
Charles's particles end up in the corresponding Bell state. To 
check whether the entanglement swapping succeeded, we as- 
sume that Alice and Charles can perform one out of two mea- 
surement x, z G {0,1} with binary outputs a, c € {0,1} 
on their system. This is sufficient, e.g., to test the CHSH 
inequality ifllll (or more precisely, for each state prepared 
by Bob, a different version of the CHSH inequality corre- 
sponding to a relabeling of the inputs and outputs). This 
scenario is characterized by the probabilities PQ(abc\xz) = 
PQ(b)PQ\i,(ac\xz), where Pq^(<ic\xz) denote the correla- 
tions seen by Alice and Charles conditioned on Bob's output 
b, and where for convenience we omitted Bob's input y since 
he is assumed to make a single, fixed measurement. 

We are interested in the robustness of Pq to the admixture 
of white noise, quantified by the maximal v € [0, 1] such that 
Pq{v) = vPq + (1 — v)Pr is bilocal, where Pr denotes the 
distribution with completely random outcomes. The quantity 
v can also be interpreted as the experimental visibility of the 
final entangled pair shared by Alice and Charles. 

If Alice and Charles use the measurement settings optimal 
for the CHSH inequality, given by xq = a x , x\ = a z , zo = 
{<j x + a z ) /y2, and z\ = (a x — a z ) /y/2, no improvement 
is obtained over the usual locality condition, i.e., we found 
that the quantum correlations become bilocal for a visibility 
v = ly/2, the same point at which they also become local. 
Using the characterization defined by and ©, we looked 
numerically for other choices of measurement settings for Al- 
ice and Charles and the best noise resistance that we found 
is v = 1/2 and is obtained for xq = zq = (a x + a z )/\/2, 
x\ = Z\ = (<r x — er z )/\/2- The corresponding correlations 
observed by the three parties are given by Pq (b) = 1/4 for 
all b, and 

!l/2 if a © c = 6 and x © z — bi 
if a © c 7^ bo and x © z = bi (6) 
1/4 otherwise 

For instance, if Alice and Charles end up in the |$ + ) state, 
corresponding to b = 00, they obtain perfectly correlated 
results if they performed the same measurements, and com- 
pletely uncorrelated results otherwise. The above correla- 
tions are local as they can be decomposed as Pq = {Pq + 
Pc)/2 where Pc and Pc are defined in term of deterministic 
strategies as P c = \ Y Ja p of 3 1 p ^Pi3 p 1 P( a ®p ){ a ®i3 ) and 

Pc = g Eaftft -P(ae/3l)(a©/3iei)-f > /3o/3i-P(a©/3o)(ae/3 ©+l)- 

They are not bilocal, however, as we now show. 

A quadratic Bell inequality for bilocality. Let us first de- 
fine, for a given probability distribution P(abc\xz), the fol- 



lowing correlation term between Alice and Charles' outputs, 
conditioned on Bob's output: 

Eb(xz) = Pb(ac\xz) — Pb(ac\xz) . (7) 

Inspired by the properties (|6]l of Pq, we introduce the fol- 
lowing combination that quantifies the high degree of correla- 
tions expected between Alice's and Charles's outcomes when 

x © z = b\ 

I = J2 P ( b ) Yl Eb ^- (8) 

b xQ)z—bi 

We quantify the deviation from the expected uncorrelated re- 
sults when x (J) z ^ bi through 

E = max max A\P(b)E b (xz)\ . (9) 

For Pq, one gets 7 = 2 and E = 0. Note that there exist bilo- 
cal correlations which also reach the value 1 = 2, for instance, 
the deterministic point defined by a = 00, (3 = 00, 7 = 00, or 
the correlations Pc introduced above. For these correlations, 
E = 4 and E = 1, respectively. The linear expression I can- 
not therefore be used as a standard Bell inequality to test bilo- 
cality. However, as we prove below, the following quadratic 
inequality 

I<l + E 2 (10) 

is satisfied by all bilocal correlations and is violated by the 
quantum point Pq, since we find 2 > 1 + 0. Intuitively, I 
correspond to a Bell inequality whose bound is not fixed for 
the entire bilocal set, but depends on how much the outputs 
are uncorrelated when x © z 7^ 61, as quantified by E. 

The noisy point Pq{v) yields I(v) = 2v and E(v) = and 
thus violate (TTOb whenever v > 1/2 (see Fig. [TJ. On the other 
hand, Pq(v) is bilocal when v < 1/2 1 1-411 - Our inequality 
thus detects optimally the resistance to noise of the point Pq. 

Proof of Ifldl. Let P be a bilocal probability distribution 
with decomposition (@}, where a = ctoax, (3 = 0of3i, 7 = 
7071. Let us define the following weights (where ao = cto © 1 
and similarly for the other indices): 

9ao«i./3o/8i.7o7i — (l a oai ,0o0i ,7o7i 9aoQi ,PoPi ,7071 

)/4, 

( ?a ai,/3o/3i,7o7i — (Saoai ,/3o/3i,7o7i ^a ai ,0o0i ,7i7o 
"^aiao,/3o/3i,7o7i ~^ 9aia ,/3o/3i ,71 7o 

The "depolarized" correlation P" = Y^a^lap-fPaPpPi 
is then also bilocal, i.e., q 1 ^ = q'^q"- Moreover, P" 
is such that I" = I and E" < E. Due to the sym- 
metries imposed through the above equations, the weights 
1a/3j depend on only 4 parameters, which we choose to be 

r = 2a =ai' S = ^7o=7i'^ = ^A)=0|aoa 1 =00,7o7i=00 m " 

u = q" a a nl nl (with obvious notations, the 

y /3 =/3i|a ai=01, 7o7l=01 v 
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FIG. 1: Two-dimensional slice of the correlation space, which con- 
tains the points Pq , Pc , Pc and Pr defined in the text [this is pre- 
cisely the slice that contains all the depolarized correlations P" in- 
troduced in the proof of J 10b ; they can indeed be written as P" — 
XPc + YPc + (1 - X - Y)Pr]. The square delimits the local 
polytope in this slice. All these local correlations can be reproduced 
in quantum theory with two independent sources. They cannot, how- 
ever, be all reproduced locally with two independent sources. The 
four portions of parabola delimit the bilocal set [the upper parabola 
is obtained from l |10| >, A similar constrained Bell-type inequality can 
be derived to obtain the lower, left, and right parabolas]. The quan- 
tum point Pq enters the bilocal region for a visibility v < 1 /2. 



weights 9„ / 3 7 being understood as probabilities 11510 . Defin- 
ing X = rs{2t - 1) and Y = (1 - r)(l - s){2u - 1), we 
find I" = 2{X + Y), E" = \X - Y\. From their defi- 
nition, X and Y are restricted to satisfy y\X\ + \J\Y\ < 
y/rs + — r)(l — s) < 1. One can easily check that un- 
der this constraint, I" < 1 + E" 2 , which implies (TTOb . 

Note that for any value of E < 1, the bound is tight, i.e., 
there exists a bilocal correlation such that I = 1 + E 2 . Take 
for instance P", with r = s = (1 + E)/2 and t = u = 1. □ 

Discussion and open questions. We have shown that if 
one makes the reasonable assumption, underlying all of mod- 
ern empirical science, that the world is composed of differ- 
ent parts that are independent for the purposes at hand, then 
nonlocality is a phenomena even more common than usually 
thought. While the standard analysis leads to the conclusion 
that the final singlet pair in an entanglement swapping ex- 
periment needs a visibility higher than v = 1 / a/2 ~ 71% 
to violate the CHSH inequality and will not violate any Bell 
inequality (with von-Neumann measurements) for visibilites 
smaller than v ~ 66% 11211 . we have shown here that under 
our assumption it already exhibits nonlocal correlations for 
visibilites as low as 50%. This simplifies the requirements for 
the demonstration of quantumness in entanglement swapping 
experiments |13 ] . 

Is this v = 50% limit a fundamental limit? It is easy to 
show that there exists a bilocal model for visibilites lower 
than 25% ifioll . But what happens in-between? Can we lower 
the visibility threshold by considering experiments with more 



inputs? Do we gain by letting Bob choose between two or 
more measurements, e.g., between two Bell state measure- 
ments in different bases? Can bilocality be violated for vis- 
ibilites lower than 33%, corresponding to a final noisy singlet 
pair that is separable? This last question is not completely 
trivial at first sight: a setting where the source Si produces a 
singlet state, 52 the separable state p = (\+z, +z) (+z, +z\ + 
\+x, — z)(+x, — z|)/2, where Bob performs a standard Bell 
measurement, Alice measures in the (a x ±cr z )/ \[2 bases and 
Charles always measures in the z basis generates correlations 
that are non-bilocal flitl . This shows that a Bell measurement 
can correlate independent systems in ways that are even more 
astonishing that one would "quantum naively" think. From 
this perspective, it would be interesting to characterize the 
class of states (including separable states) that are non-bilocal 
when correlated through Bell measurements. 

The present work raise many other questions. In particular, 
the condition (0) can be straightforwardly extended to models 
with n independent sources. How do such n-local models dif- 
fer from bilocal ones? How does their tolerance to noise scale 
with nl Finally, it would be interesting to explore the implica- 
tions of our approach and findings in the context of quantum 
information protocols based on non-locality, and in particular 
protocols that use at their heart measurements on independent 
systems, such as quantum repeaters and measurement based 
computation. 
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Photonics and the European ERC-AG QORE. 
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r,s,t,u £ [0,1]. The weights <j^a~ are then fully defined by 
n" — -fi - ■ a" — -S I - • 

TOOK! 2 Ua 0, a l^ 2 Ua 0,°'l ' 'iJoJl 2 "70>7l~ 2 u 70.71' 

Ift*-, = |< 5 Qoe7o i /3o + i T i ^o©7o,/3o if Q o = Q i and 70 = 71; 
9/3| a7 = f<^oe7o,A)©/3i + ^ ii ^oe7o,/3o©/3iei if a <) / Q l 
and 70/71; q'p\ an = j otherwise. 
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CHSH inequality performed on the singlet Si, where Charles's 
outcomes specify in which basis (x or z) Bob's particles has 
been measured. Such an interpretation makes the similarity be- 
tween our independence assumption and the "free choice" as- 
sumption of standard Bell experiments more explicit JToll . 



